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Abstract 

After n random polarizations on S d , the expected symmetric difference of a Borel set from a 
polar cap is bounded by C^n , where the constant depends on the dimension [1]. We show 



o 

as: 
< 
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psj , here that this power law is best possible and that necessarily C,i > d. 

Ph 1 Introduction 

Let A be a subset of the rf-dimensional sphere E> d (viewed as the unit sphere in and let a : x M- x 

be a reflection at a great circle that does not pass through the north pole. The polarization of A with 
respect to a is defined by 

> 

SA f x e A or x e A, if S(x,0) < 6{x,0) 
>T~j ■ 1 .t G A and x G A , if 5(x, O) > Six, O) 

in 



Here, y) denotes the geodesic distance on E> d given by the angle enclosed between x and y, and O 
denotes the north pole. Since reflections preserve the uniform probability measure m(-) on the sphere, 
so do polarizations, and 

m(SAnSB)-m(Af]B)= / I A \b(x) I b \a(x) dx > . (1) 

Js d 

We parametrize the reflections on E> d by w G = § d / ±, setting 

&u(x) = x — 2(u ■ x)u , 

and we denote the corresponding polarization by S u . A random polarization Su is polarization in 
the direction of a uniformly distributed random variable U on f2. We consider sequences of random 
polarizations Su x „,u n = Su n o- ■ -oS^, where the {t/j}j>i are independent. Van Schaftingen has shown 
that almost surely, for every Borel set A the sequence S^.^UnA converges to the polar cap A* of the 
same volume [5, Theorem 3.13]. The convergence occurs in symmetric difference if A is measurable, 
and in Hausdorff metric if A is compact [1 , Corollary 4]. 
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Subject of this note is the rate of convergence. In prior work, we have shown that under a similar 
sequence of random polarizations on ~§t d , 



There, A is a Borel measurable subset of the unit ball, m is Lebesgue measure (normalized so that 
the unit ball has measure one), and Cd = d2 d+1 03 Proposition 4.1]. This rate of convergence is 
much slower than what is known for other symmetrizations. Klartag has proved that a sequence of 
3d carefully chosen Steiner symmetrizations in R. d followed by a random sequence where each step 
consists of d orthogonal Steiner symmetrizations converges faster than every polynomial. The leading 
constant depends only on the dimension and grows at most polynomially [4, Theorem 1.5]. Although 
Klartag's result applies only to convex bodies, it raises the question whether the power law in Eq. © 
can be improved. For random polarizations, the answer is negative: 

Proposition. For random polarizations of a Borel set A C E> d , Eq. © holds with Cd = 2 d . The rT x 
power law is optimal, and the sharp constant satisfies Cd > d. 

The proof of the proposition has two parts. Eq. © and the upper bound on Cd are obtained by 
simply adjusting Proposition 4. 1 of [[II to the sphere. For the lower bound on Cd and to prove the opti- 
mality of the power law, we analyze how spherical caps move under polarization. If A is a hemisphere, 
we compare the difference of Su x ...v n A from A* with with the order statistics of the uniform distribu- 
tion, and relate the limiting distribution of n ■ m(Su 1 .„u n A A A) to a Gamma distribution. We work 
on the sphere rather than on M. d , because the additional symmetry simplifies the calculations. It will 
be clear from the proofs that similar bounds hold for the polarization of balls in M. d . Other questions 
remain open: How quickly do the sharp constants grow with the dimension? What is the impact of the 
distribution ofU? Can one speed up the convergence by adapting the sequence to A? 

Acknowledgments. I wish to thank Vitali Milman for his comments on [1] that prompted this work, 
and Marc Fortier for patience and helpful discussions. The research was supported in part by NSERC 
through Discovery Grant No. 311685-10. 

2 The upper bound on the sharp constant 

In this section, we show that Eq. © holds on the sphere with Cd = 2 d . For a single random polarization 
we have by the identity © and Fubini's theorem, 



We compute the probability under the integral as an average over the hemisphere where u ■ x > 0, and 
change variables to z = a u (x). The inverse of the map u M- z and its Jacobian on the tangent space of 
§ d are given by 




(2) 




u(z] x) 



x — z 



3 {z-x) = {2 




x — z 
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where \x — z\ = 2 sin <5 ^' 2 - ) is the Euclidean distance between x and z in M. d+1 . We obtain 

m(AAA*)-E[m(SuAAA*)] = 2/ / |x - zl"^ 1 ) dm{z)dm{x) 

Ja*\a Ja\a* 

> 2- d (m(A A A*)) 2 . 

For a random sequence Sux...u n , we take expectations again and apply Jensen's inequality to see that 



E 



m(S Ul ...u n - 1 AAA*)-m(S Ul ...u n AAA*j\ > 2~ d (E [m(S Ul ... Un _ 1 A A A 



It follows that z n = 2~ d E[m(S Ul ... Un A A A*)] satisfies z~ l > z n \ + 1, proving Eq. © with constant 
C d = 2 d . □ 



3 Random compressions 

Let A be a spherical cap centered at a point a. Polarization with respect to a reflection a : x h->- x 
transforms A into the spherical cap of the same volume centered at r(a), where 



T(X) 



x, 6(x,0)<6(x,0) 
x , otherwise . 



We will refer to r as the compression associated with a. The compression associated with a random 
reflection o\j will be denoted by tjj. The following lemma describes the distribution of the distance of 
Ty{x) from the north pole. 

Lemma. If U is uniformly distributed on Q, then for every point x G E> d with S(x,0) = £ 

P(5(r u (x),0)>P) = I^Jl-- /7 COS ^~ C ° S ^ ) (d ' 1)/2 ^| , /3e[0,7r]. (3) 

Proof. By definition of the compression, 

P(8(tu(x), O) > 0) = I i>p P(5(au(x), O) > P) . 

For £ < /3, there is nothing more to show. For £ > (3, we set t = cos /3 and calculate the spherical 
average as an expectation with respect to the standard normal probability measure on M. d+1 , see B2l 
Exercise 63 on p. 80]. We use the coordinate system u = (r cos 0, r sin 6, u) Gixlx R d_1 ,which 
we rotate into a position where x = (cos |, sin |, 0) and O = (cos |, — sin |, 0) . Then 

(u ■ x){u ■ O) = — (cos 2$ + cos£) , 
and 5(cr u (x), O) < p if and only if — r 2 (cos26 ) + cos£) > (\u\ 2 + r 2 )(t — cos£). This results in 

P(S(tu(x),0) < p) = / I au{x y >tdm(u) 

= 7T / / / ^-r 2 (cos2e+t)>|«| 2 (t-cosO 2re" r dr d>y(u) d6 , 
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where 7 is the standard normal probability measure on M rf 1 . We integrate explicitly over r and 

t— cos ^ N 
t+cos 20 / 



evaluate the remaining Gaussian integral by rescaling v = (l — / +c c ° s 2 g ) ^ 



P(5(ra(x),0) </3) = ^ f J cos20+t < o ( / e^ 1 " 1 ' dj{u)\ d6 

1 r /-cos 20 - cos/3\(^-i)/2 

— — / -'cos26>+t<o( ^ r J 

Z7T 7-tt V — COS 2u — COS £ / 

The claim follows after restricting the integral to a half-period and changing variables 29 — > n — 6. □ 
For d = 1, the reflected point ajj{x) is uniformly distributed on S 1 , and Eq. © reduces to 

P(8(ru(x), 0)>(3) = Is (x ,o)>(s (1 - f ) , /? G [0, vr] . 

As ti increases, ^(a;) concentrates in a ball of radius comparable to d^ 1 / 2 about x, its distance from 
the north pole concentrates in an interval of length comparable to d^ 1 about £, and the integral in 
Eq. © goes to zero. For all d > 1 and < /3 < £ < 7r, we have the bound 



1 ^/ cosg-cos/? y^)/2^ < g / l- cos/3 x(rf-i)/2_ | 2 sin f /sin 

7r 7 V cos 6* — cos £/ ~~ 7r V 1 — cos £ / 7T \ sin 



If {t/j}j>i is a sequence of independent uniformly distributed random variables in fi, it follows that 
the Euclidean distance Y n = \ru 1 ...u n (x) — 0\ satisfies the recursion 

P(Y n+1 >r 1 \Y n )> I Yn>v {l - jly-)^ 1 } , [0, t\ (4) 

with initial value Y = \x — 0\ = 2 sin | and with £ = tx — 0(£ 2 ). 



4 The lower bound on the sharp constant 

Let A be the hemisphere centered at a point a ^ O, and set a = 5(a,0). We claim that 

liminf nElmiSu, Un A A A*)] > (1 - 0(a 2 )) d . (5) 

Taking a — > 0, we see that the sharp constant satisfies C d > d, completing the proof of the proposition. 

To prove the claim, consider a sequence of random points {Vi}i>i that are distributed independently 
and uniformly on an interval [0, £], and let Y n be the <i-th lowest point among Vi, . . . , V n+d . The random 
variable Y n is called the d-th order statistic of V\ . . . , V n+d . The sequence {Y n } n > solves Eq. © 
with equality, because conditioned on Y n = y, the d— 1 points among V\, . . . V n+d to the left of y are 
independent and uniformly distributed on [0, y), and V n+d +i is independent and uniformly distributed 
on [0, £}. The joint distribution of the order statistics can be written explicitly in terms of binomial 
random variables B(n,p), see (31 Exercises 21-25 on p. 142]. We have 

P(Y n >ri\Y = y) = I y>v £ P(B{d-l, J) = j) ■ P(B{n, ?) = *), 

j+k<d 
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where the first factor in the sum accounts for the points among Vi,...Vd that fall to the left of r\, while 
the second factor accounts for such points among Vd+i, . . . V n+ d- By Stirling's formula, 

P{nY n > v \Y = y) P (T(d) > f) ( n oo) 

for each y G (0, £], where r(d) is a Gamma random variable that describes the d-th point in a Poisson 
process of intensity one Exercise 24 (b) on p. 142]. In particular, E[F n | F = 2/] — > 2d. 

The center of S^., .u n A is given by Tu x ...u n (o). We have shown in Section 3 that F n = Ir^...^ (a) — 
0\ satisfies Eq. ©. Since the right hand side of this recursion increases with Y n and the geodesic 
distance on the sphere exceeds the Euclidean distance, 

P{5{Tu,...u n {a), 0)> V )> P(Y n > v \Y = a) 

for all n > with £ = tt — 0(a 2 ). For the mean, this implies that 

liminfnE[5(7v 1 ...ir n («),0)] > (n-0(a 2 ))d. 

n— >oo 

Eq. © follows because the symmetric difference between two hemispheres is just the distance of their 
centers, expressed as a fraction of tt. □ 

Remark. A slightly more careful analysis of Eq. © shows that for a ^ O, 

nS(T Uu .. Un (a),0) -> irT(d) (n -> oo) 
in distribution, and hence lim n F,[m(Su 1 ...u n A A A*)} = d for the hemisphere A centered at a. 
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